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Electromagnetic Structure of Light Baryons in Lattice QCD 
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A method in which electromagnetic properties of hadrons are studied by direct simulation of dynamical photon 
effects is applied to the extraction of the isomultiplet structure of the octet baryons. Using 187 configurations at 
f3 — 5.7 with Wilson action, and up and down quark masses determined from the meson spectrum, the nucleon 
splitting is found to be 1.55(±0.56 stat) MeV; the hyperon splittings are found to be E° - E + = 2.47 ± 0.39, 
£~ — S° = 4.63 ± 0.36, E~ — H° = 5.68 ± 0.24 MeV. Estimated systematic corrections arising from finite volume 
and the quenched approximation are included in these results. 
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1. Introduction 

The fact that the intrinsic up-down quark mass 
difference is comparable to typical hadronic elec- 
tromagnetic energy shifts has made it difficult to 
reliably calculate the octet baryon isomultiplet 
splittings in the absence of a systematic treat- 
ment of virtual photon effects combined with non- 
perturbative QCD contributions. The problem 
is particularly acute in the case of the proton- 
neutron mass difference, where tadpole (quark 
mass) effects almost completely cancel the pho- 
ton cloud contributions (see [Q for an early re- 
view exhibiting the massive confusion prevailing 
in pre-QCD days). In this talk, we extend a 
method recently used Q to extract electromag- 
netic contributions to pseudoscalar masses to the 
octet baryon spectrum. The basic idea is to prop- 
agate the quarks through a U(3) = SU(3) x U (1) 
field including both dynamical gluon and pho- 
ton effects. The calculation of baryon isomultiplet 
splittings can then be viewed as a two-step pro- 
cess: 

(a) First, the bare quark masses are determined 
from an analysis ||] of the pseudoscalar me- 
son spectrum (including electromagnetic contri- 
butions); 

(b) Secondly, the octet baryon spectrum is com- 
puted (again including dynamical photon effects) 
and extrapolated to physical values of quark mass 



(as determined in step (a)) and electric charge. 

2. Extraction and Fitting of Baryon Spec- 
trum 

The strategy of the calculation is as follows: 
quark propagators are generated in the pres- 
ence of Coulomb gauge background SU(3)xU(l) 
fields. 187 gauge configurations, separated by 
1000 Monte Carlo sweeps, were generated at fi = 
5.7 on a 12 3 x 24 lattice. Quark propagators are 
calculated for 4 electric charges and 3 light quark 
mass values, and with either a local or smeared 
source (see for details). From the resulting 12 
quark propagators, 936 independent octet baryon 
three-quark combinations can be formed. 

In quenched QCD it is known || that baryon 
masses are described by a function of the bare 
quark masses involving nonanalytic m 3 / 2 (as well 
as linear) terms, and terms involving logarithms 
of the quark mass arising from the same hairpin 
diagrams familiar in the quenched meson spec- 
trum [HjU . The latter terms now appear to be ex- 
tremely small numerically ||[7| : we neglect them 
throughout. However, we do include terms of or- 
der (quarkmass) 3 ^ 2 . Thus a general octet baryon 
mass is written 

tub = A{e q i,e q 2,e q z) + ^^m, ql Bi{e q \,e q 2,e q s,) 
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(m gi +m qj ) 3/2 C ij (e q i,e q2 ,e q3 ) (1) 



where e 9 i , e g 2 , e g 3 are the three quark charges, 
and m q i , m q 2 , m g 3 are the three bare quark 
masses, defined in terms of the Wilson hopping 
parameter by (k^ 1 — n~ 1 )/2a. (Here a is the lat- 
tice spacing.) Each of the coefficients A,Bi,Cij 
in (|l|) is then expanded in powers of the quark 
charges e q i, e q 2, e g 3, with terms up to fourth order 
for A, second order for Bi, and with no charge de- 
pendence assumed for the nonanalytic CV,- terms. 
(|l|) turns out to have 30 parameters once all sym- 
metries arc exploited. 

We have varied the baryon mass window (for 
each choice of Euclidean time window used to ex- 
tract a mass from smeared-local correlators) until 
the x 2 /dof of the fit to (|l|) was minimized. For ex- 
ample, using a Euclidean time window from t =5 
to t =8, the mass window (lattice units) from 
1.20 to 1.26 was found to contain 74 baryons. 
Determining the 30 parameters in ([!]) by fitting 
this set of masses gave a x 2 /dof of 1.33 . By 
contrast, using the mass window from 1.15 to 
1.20 (122 baryons), the chi-square fit minimizes 
at x 2 /dof=2.16. For each choice of Euclidean 
time window, we have performed the fit to ([!]) 
using a baryon mass window which optimizes the 
X 2 /dof - the overall optimal fit was found for the 
window t = 5 to 8. One then determines the 
mass of any given octet baryon by extrapolating 
to physical values of quark mass and charge. The 
propagators for different electric charge are highly 
correlated, so it is not surprising that the statisti- 
cal error on the center of gravity of isomultiplets 
is considerably larger than the error on multiplet 
splittings. Note, in connection with the raw lat- 
tice results quoted in Table 1, the following : 
(1) The quark mass parameters and lattice scale 
assumed in generating masses for each of the fit- 
ting window choices are shown in the second row. 
The up and down quark masses are those ob- 
tained from the pseudoscalar spectrum (these will 
depend on the scale). The strange quark mass is 
known to fall at a higher value when determined 
from the baryon spectrum (due to discretization 
and quenched errors), so we have chosen to fix 
it using the center of gravity of the 2 hyperon, 



Table 1 

Raw lattice results (McV) for baryon octet 
(/3=5.7, 12 3 x24, 187 configurations) 



Baryon State 


Window 5-8 (x 2 /dof=1.33) 


Parameters 


m u ,d,s,a L =3. 57, 7. 10,155,1370 


N 


935.92±42.4 


P 


933.07 ±42.9 


N-P 


2.83 ± 0.56 


E+ 


1171.6 ± 25.6 


E° 


1175.1 ± 25.3 


sr 


1179.1 ± 25.0 


E°-E+ 


3.43 ± 0.39 


E" -E° 


4.04 ± 0.36 


E+ + E" - 2E° 


0.61 ± 0.19 




1312.9 ± 14.5 


wo 


1308.2 ± 14.6 




4.72 ± 0.24 


" A " 


1098 ± 52 



which has the smallest statistical errors in our 
analysis. The center of gravity of the E multiplet 
and the A mass are then predictions of the anal- 
ysis. 

(2) The lattice scale has been fixed in each case 
by requiring the nucleon center of gravity to sit 
at (roughly) the physical value. 



3. Finite Volume and Quenched Correc- 
tions 

With massless physical degrees of freedom we 
expect finite volume effects which fall as inverse 
powers of the lattice size. These corrections can 
be studied directly on the lattice by repeating the 
calculations on lattices of varying physical vol- 
ume. Here we estimate them by using the known 
dominance of the Born contribution to the disper- 
sive evaluation of the Cottingham formula. Sin- 
gle photon exchange can be written as a sum of 
an electric and magnetic contribution to hadronic 
self-energies- the electric term takes the form 

c 9 1 \- G E (q) 2 r 2 

Sm el = 2,am-^- 1 ^{- r ^- 2 
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Table 2 



Final results for baryon octet splittings ((3=5.7, 12 3 x24, 187 configurations) 



Level Splitting 


Raw Lattice 


Finite Volume 


Meson Cloud 


Total Lattice 


Physical 


N - P 


2.83 ± 0.56 


-0.75 


-0.53 


1.55 ± 0.56 


1.293 




3.43 ± 0.39 


-0.80 


-0.16 


2.47 ± 0.39 


3.18 ± 0.1 




4.04 ± 0.36 


+0.86 


-0.27 


4.63 ± 0.36 


4.88 ± 0.1 


£+ + IT - 2S° 


0.61 ± 0.19 


+1.66 


-0.11 


2.16 ± 0.19 


1.70 ± 0.15 




4.72 ± 0.24 


+0.86 


+0.10 


5.68 ± 0.24 


6.4 ± 0.6 




-1)} (2) 



where the momentum vectors q are the discrctizcd 
bosonic photon momenta for the finite LxLxL lat- 
tice. Using (^|) one can estimate the finite vol- 
ume corrections to baryon masses on our L=12 
lattice- they are indicated in column 3 of Table 
2, together with our final estimate (including the 
finite volume correction as well as quenched er- 
ror estimate- see below) for the baryon mass in 
column 5. 

Processes in which mesons are emitted and re- 
absorbed from a baryon include graphs with in- 
ternal quark loops and are known || to result 
in a small but nonnegligible shift in isospin split- 
tings. For example, in the static limit where the 
nuclcon mass is infinite, the pion cloud decreases 
the neutron-proton splitting by an amount (in 
the infinite volume limit) 0.43AMo, where AMo 
is the nucleon splitting in the absence of a vir- 
tual pion cloud (a fully relativistic evaluation 
gives 0.41AMq). We shall use a static approxi- 
mation but include the effects of all octet pseu- 
doscalar mesons (assuming SU(3) symmetry with 
a d : (f + d) ratio of 0.62). Discretizing the second 
order shift formula (see |) on a LxLxL lattice, 
one may estimate the meson cloud shift for the 
particular lattices used. Since the meson cloud 
shift includes contributions from quenched non- 
planar graphs in the cases where the emitted me- 
son only contains valence quarks of the external 
baryon, these estimates are only a rough indica- 
tion of the magnitude and sign (probably, an over- 
estimate), of the quenched correction. Setting 
L=12 and using a lattice scale a~ 1 = 1370 MeV, 
together with the quenched masses from column 2 
of Table 1, we obtain the meson cloud shifts given 



in column 4 of Table 2. The lattice results, cor- 
rected for finite volume and meson cloud effects, 
are given in column 5, and the physical values in 
column 6. 

The results in Table 2 (which must still be cor- 
rected for finite lattice spacing effects) suggest 
that this first evaluation, on a fairly coarse lat- 
tice, already reproduces - almost quantitatively- 
the isomultiplet pattern of the octet baryons. Of 
course, the extremely delicate level of baryon fine 
structure being considered here requires a de- 
tailed study of all systematic effects, with im- 
proved statistics on larger lattices. An upcoming 
run will work with lattices of varying physical vol- 
ume, using improved (to 0(a)) action to minimize 
lattice discretization errors- an important check 
given the known strong a-depcndcncc in off-shell 
defined continuum quark mass parameters ||. 

We thank George Hockney for continuing con- 
tributions to our effort. 
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